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ABSTRACT 



pL| I Bernardeau et al. (1997), using perturbation theory, showed that the skew- 

0\ [ ness of the large-scale lensing-convergence, or projected mass density, could be 

used to constrain Q^n, the matter content of the universe. On the other hand, 
deep weak-lensing field surveys in the near future will likely measure the conver- 
if^ • gence on small angular scales ( "z 10 arcmin.), where the signal will be dominated 

!>■ 
(N 

CN ' scale convergence skewness, making use of a prescription for the highly nonlinear 



three-point function developed by Scoccimarro and Frieman (1998). This method 



0> ■ gives predictions that agree well with existing results from ray-tracing N-body 

1^ • simulations, but is significantly faster, allowing the exploration of a large num- 

I ' ber of models. We demonstrate that the small-scale convergence skewness is 

;_( ■ insensitive to the shape and normalization of the primordial (CDM-type) power 

^ ■ spectrum, making it dependent almost entirely on the cosmological energy con- 

tents, through their influence on the global geometrical distances and fluctuation 
1^. growth rate. Moreover, nonlinear clustering appears to enhance the differences 

j_j ■ between predictions of the convergence skewness for a range of models. Hence, 

in addition to constraining fi^, the small-scale convergence skewness from future 
deep several-degree-wide surveys can be used to differentiate between curvature 
dominated and cosmological constant (A) dominated models, as well as to con- 
strain the equation of state of a quintessence component, thereby distinguishing 
A from quintessence as well. Finally, our method can be easily generalized to 
other measures such as aperture mass statistics. 



Subject headings: cosmology: theory — gravitational lensing — large-scale struc- 
ture of universe 
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1. Introduction 



The correlated shear of images of distant galaxies provides a promising way to probe 
the intervening large scale structure of the universe (e.g. [B landlord et al. 199 1| ; pVliralda-] 
Escude 1991 ; Kaiser 1992 ; see also |Jain et al. 1999| [JSW hereafter] and ref. therein). The 



convergence can be constructed from a shear map ([Kaiser fc Squires 19931) , which can be 
interpreted as a form of projected mass density ( [Kaiser 1992| our notation follows that of 
[Jain fc Sel.iak 1997| ): 

K{e)= rdxw{x)S{r{x)0,x) (1) 

Jo 

where 6 is the mass overdensity as a function of spatial position (and implicitly as a function 

of time as well, where the time and space coordinates fall on the photon null geodesic), 

6 is the angular position on the sky, x is the comoving distance along the line of sight, 

r{x) is the angular diameter distance, and w{x) is a weight function which depends on a 

combination of global geometrical distances and is proportional to the total matter density 

of the universe. The coordinates x = , Xs denote respectively the positions of the observer 

and the sources or background galaxies. The so-called Born approximation has been assumed 

(see [Bernardeau et al. 1997[ for a discussion). 



It is clear that k is a valuable quantity for cosmology. One can derive important con- 
straints on the cosmological density parameters fi's through the dependence of n on the 
global geometrical distances and the evolution of 6 on the line cone. A commonly used 
statistic is its second moment (k^), or the two-point correlation function {K{d)n{d')). How- 
ever, it is clear that the second moment depends on the mass power spectrum as well (e.g. 
Jain k SeljakTW^ ; [Kaiser W^ . 



Bernardeau et al. ([1997[) , using perturbation theory, showed that this degeneracy could 
be broken by using the convergence skewness 5*3 = {k^) / {k^)"^ . It is customary to consider 
S's as a function of angular scale On, assuming n is first smoothed on scale Or. However, 
future weak lensing surveys are likely to yield measurements of 5*3 first on small angular 
scales. Or < 10', both because the small-scale shear signal is stronger and also because of 
larger sampling fluctuations on large scales (see e.g. [Van Waerbeke et al. 1999| , JSW) . For 



sources at a redshift of z = 1, the peak contributions to the lensing signal will come from 
z ~ 0.5, which for Or < 10' translates into a comoving length scale of less than a few Mpc 
that is generally comparable to or smaller than the nonlinear scale (where the rms density 
fluctuation is of order unity). It is therefore expected that the perturbative treatment of 
Bernardeau et al. ( [1997[ ) would not hold for these angular scales of interest. This has in fact 
been explicitly demonstrated by JSW using the technique of ray-tracing N-body simulations 
(see their Fig. 18; see also [Couchman et al. 1998| ). Unfortunately, the prediction of skewness 
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from N-body simulations can become prohibitively expensive, if one is interested in exploring 
a large number of cosmological models. 

There is therefore a need for alternative methods to predict accurately and efficiently 
the small-scale skewness. It is our aim here to develop such a method. We will test it by 
comparing with existing results from simulations, and show that the small-scale skewness is 
a sensitive probe of fi's. We will then apply it to cosmological models that have not been 
considered before in the context of weak lensing. In particular, we will predict the small 
angular-scale skewness for a quintessence model, where quintessence is a component of the 
cosmological fluid that has negative pressure (e.g. [Peebles fc Ratra 1988| ; [brieman et al. 1995 



Coble et al. 19^; frurner fc White 1997|; [l^'erreira fc Joyce 1998|; ^^aldwell et aLT^gSQ. Such 



models, which include the cosmological constant dominated models as a limiting case, are 
currently in favor in part because of recent Type la supernova measurements ( [liiess 19^ 



Garnavich 1998; Perlmutter 1998). We will demonstrate that the convergence skewness can 



provide interesting constraints on them. 



2. The Convergence Skewness 

Let us first give the expressions for the cosmology dependent geometrical quantities that 
appear in eq. (|^). The comoving distance along the line of sight x is given by (Peebles 1993|) 



X{z) = cH^' r dz'[nU^ + zf + fifc(l + z'y + fi,(l + ^')3(i+-.)]-i/2 (2) 

Jo 

where z is the redshift of interest, Hq is the Hubble parameter today, c is the speed 
of light, and the fi's denote the fractions of the critical energy density today in vari- 
ous components: Qm for pressureless matter or dust, Qk for spatial curvature and fig 
for quintessence or a fiuid with negative pressure (its pressure p is related to its den- 
sity p hj p = Wqp, where Wq < 0), with the cosmological constant A as a limiting case 
{wq = —1). The fi's sum to unity. The angular-diameter distance r{x) is given by 
r(x) = K"^^'^ sin K^^'^Xi Xj (—K)^^^"^ smh(—Ky^'^x for closed, fiat and open models respec- 
tively, and K = {Qm — ^)c^^Hq. In other words, the metric is given by ds^ = —c^dt'^ + 
O'i'tyidx^ + ^(x)^'^^^)) where ait) = 1/(1 + z) is the expansion scale factor as a function of 
proper time t. The line-of-sight projection of 6 in eq. (|l]) is weighed by the function w{x)'- 

where Xs is the comoving radial position of the sources. Note that here, as in the rest of 
the paper, we assume all sources are at the same redshift. Eq. (|1]) and (|^) can be easily 
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generalized to the case of multiple source-redshifts by integrating over contributions from 
different Xs's- Statistical measures of k for sources distributed in a realistic fashion can 
usually be approximated by having all sources at the same mean redshift (e.g. JSW). 



The convergence skewness is defined by 

s,{eR) ^ p^^ , ^9n^ f <e')Wg,{e - 0')<fe' (4) 



where Wg^ is a smoothing kernel of radius Or (in this paper, we will use a top-hat). The 
utility of S'3 derives from, crudely speaking, the fact that its analogue for the mass overdensity 
(((5^)/(5^)) is quite insensitive to details of the power spectrum, especially on small scales. 
Hence, as we will see, S'3 is almost purely determined by the cosmological energy contents. 



Combining eq. (|I|) and (^, it can be shown that 

S3 = K^/iK^f with (5) 

K, ^ (27,f j^' dx^^ j dH,dH,B{l,/r{x),e-2/r{x),e-^lr{x)) 

w{heR)w{^2eR)w{i^eR) 

K, ^ 27^ j^' dx^^ j dHP{^/r{x))W{mRf 

where W{x) is the Fourier transform of the two-dimensional top-hat, W{x) = 2Ji{x)/x with 
Ji being the first-order Bessel function. The ts represent the Fourier coordinates in angular 
space, in other words we are taking the small-angle approximation where spherical har- 
monics can be replaced by plane waves. The combination £1 + €2 + -^3 vanishes. The three- 
dimensional mass power spectrum and bispectrum are respectively P and B. Our convention 
is: ^2(|r|) = J d^kP{k)e-'^'' and 6(ri, jral) = J d^kid^k2B{ki,k2,-'k:L - k2)e-*^^'"^-*^^-''^ 
where the k's denote the Fourier coordinates in three-dimensional space, and ^2 and ^3 are the 
two- and three-point correlation functions respectively. (For readers who are used to putting 
(2-7?)^ under d^k: simply replace all relevant expressions in this paper by P — >■ P/(27r)^ and 
B -^ B/{27t)^.) Note how in eq. (^ the projection forces the k's to lie in the plane of the 
sky. The reader is referred to Kaiser (|1992| ), Bernardeau et al. ( |1997[ ) and Jain & Seljak 
( |1997| ) for detailed derivations. 

To compute S3 on small angular scales, we need to understand the nonlinear evolution 
of P and B. The nonlinear behavior of P can be described by a scaling ansatz introduced 
by Hamilton et al. (|1991| ), which was later extended by Jain et al. ( |1995| ) and Peacock 
& Dodds ( [1994^ |1996| ). We will employ the latest version set out in the latter. Jain & 



Seljak ( |1997| ) have considered the two-point version of K2 using this ansatz. Essentially, the 
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ansatz consists of postulating that Airk^Plk) = /[47r/c£P/;,(fc/^)] where / is some universal 
function, and PLiki) is the linear power spectrum at the rescaled wave- number defined 
by ki = [\ + Ai^k^ P{k)]~^/^k. The cosmological dependence comes in through the linear 
fluctuation growth rate Pl oc [g{z)/{l + z)Y (fitting formula from UarroU et al. 199"^ ): 



g{z) = ^-n^{z)[n^{zYl' - n^{z) + (1 + n^{z)/2){l + fiA(^)/70)]-^ (6) 

where Vtm{z) = fir„(l + zf/[nm{l + zf + l]fe(l + z)^ + ^a] and n^iz) = f]A/[fi^(l + zf + 
r2fc(l + z^ + I^a], and fi's without explicit z dependence denote their values today. For 
quintessence models with Wg ^ —1, we integrate numerically the equation for the linear 
growth rate, and substitute this in the corresponding expressions given by Peacock & Dodds 
( |1996|) (see [Wang fc Steinhardt 1998| for a useful fitting formula). 

For the bispectrum, it has been conjectured for some time that the following scaling 
approximately holds in the highly nonlinear regime (e.g. pavis fc Peebles 198(J| ; [Peebles 198U ; 
[Fry 19841; [Hamilton 19"88|) : 

5(ki, k2, ka) = Q3iPik^)P{k,) + P{k^)P{k,) + P{k;)P{k,)) (7) 

where the three k's form a closed triangle, and Q3 is a weak function of scale but independent 
of the triangle configuration. The above nonlinear hierarchical form (sometimes called the 
hierarchical ansatz) is also observed in N-body simulations (see pcoccimarro et al. 1998[ and 



ref. therein). The problem, however, was that, there has been for a long time no way 
to predict the amplitude of Q3, other than by examining N-body simulations on a case 



by case basis. Recently, Scoccimarro fc Frieman ( |1998|) introduced a method they named 



hyperextended perturbation theory which allows one to calculate Q3 analytically: 

g3(n) = [4-2"]/[l + 2"+i] (8) 

where n is the linear power spectral index at the scale of interest (fci + k2 + k^)/?). The above 
expression implies that Qs is insensitive to cosmology, except through n (see ^coccimarro 
et al. 19981 ). 



Hence, combining eq. (0), (H) and the nonlinear evolution of P given by Peacock fc 
Dodds ([1996[) , together with eq. (H), completely specifies 5*3 for any given primordial power 



spectrum and cosmology. To ease the computation, we find that the following approximation 
for i^3 agrees with the exact integration to within a few percent for the models considered 
in this paper: 



<-Xs w{xf 



K, ~ 3{2rrf £' dx^^ [J d'i^^P{i/r{x))W{eer,, 



(9) 
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where Q3 is evaluated at an n corresponding to the scale C-/r{x). This approximation works 
in part because Qs varies slowly with scale, on the relevant small scales. 

In Fig. |l|a, we show a comparison of the skewness computed as described above with the 
skewness obtained from ray-tracing N-body simulations (JSW), for sources at z = 1. The 
error-bars shown are estimated from the dispersions between 5-10 ray-tracing realizations. 
Three models are shown (see Table p. They are all normalized to match the cluster abun- 
dance today. The agreement is good, to better than 10% for 0/j ~ 1' — 5', and it remains 
reasonable at larger angular scales, although its exact level is somewhat uncertain because 
of the large dispersions of the N-body results. The agreement here is to be contrasted with 
the as much as 30% discrepancy for the perturbation theory predictions, shown as points on 
the far left of Fig. |T|a. In particular, perturbation theory brings S^ for OCDM and LCDM 
much closer than what it should be (the good agreement of the perturbative 6*3 with the 
actual value for LCDM seems to be a coincidence). It is interesting how nonlinear clustering 
makes it easier to tell them apart. There seems to be a complicated interplay of projection 
and nonlinear clustering (e.g. |Gaztafiaga fc Bernardeau 1998|) . Our OCDM predictions seem 



to be systematically a little higher than the N-body results, but it should be kept in mind 
that measurements at different scales are correlated and that a measurement-bias due to a 
division of estimators might be present ( [Hui fc Gaztanaga 1999| ). 



The accuracy of our method is actually somewhat surprising because of the inherent 
approximate nature of the prescription for Qs (eq. |§|) and of the scaling ansatz for the 
power spectrum evolution. Moreover, the weakly nonlinear fluctuations, which do not obey 
the hierarchical form with a configuration independent Q3 as in eq. (|^, must contribute 
at some level to the relevant integrals for 6*3. To check this, we perform an alternative 
integration for K2 in eq. (^ by including only 'nonlinear' modes: taking the lower limit of 
integration to be (i^\ instead of 0, where £ni satisfies 47r(£ni/r(x))^P(£ni/r(x)) = 1. Let us 
call the resulting integral K2^, and define /S.K2 — 1-^2^ ~ K2\/K2. We find that /S.K2 is very 
similar for all 3 models above, and is about 10 — 30% at Or r^ V — 5', reaching about 45% 
at Or ~ 10'. We therefore propose the following self-consistency check: A^^^ should be less 
than about 30% for our method to yield reliable estimates of 5*3. 

We show in Fig. 0b our prediction of 5*3 for the same three models, but the points with 
error-bars now represent measurements from simulated surveys, of a size 3° x 3°, with 2 x 10^ 
z = 1 galaxies per square degree whose intrinsic ellipticities are Gaussian distributed with 
an rms of 0.4 for each component (attainable with multiple several-hour-long exposures on 
a 4-meter class telescope using large CCDs; |Van Waerbeke et al. 1999| , JSW). It is clear that 
such a survey can separate these 3 models very nicely. Note, however, systematic errors have 
not been taken into account. 
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Perhaps more interestingly, we show in the same figure our prediction of S^ for a cluster- 
normahzed ( |Wang fc Steinhardt 1998| ) quintessence model (qCDM, see Table p. The equa- 
tion of state {wq = —0.5) is motivated by certain models of dynamical supersymmetry 
breaking (e.g. Binetruy 1999| ). Fluctuations in the quintessence component have been ig- 



nored, which is probably a good approximation on the small scales of interest ([Turner fc 



White 19^ ; [Caldwell et al. 1998|) . Wang et al. ( [1999[ ) have argued that current observations 
cannot tell apart qCDM models with —l<Wq^ — 0.4 from LCDM models, for Q^ ~ 0.3. 
Future microwave background experiments would be able to provide better constraints, but 
there exist significant degeneracies, especially if Hq is allowed to vary ( [Huey et al. 1998[ ). 
The skewness has the advantage that it is independent of Hq. Fig. |l]b shows that the small 
angular-scale convergence skewness (especially at 6f( ~ 1' — 5') provides a promising way 
to disentangle qCDM and LCDM models: fixing Q^, S^ varies smoothly from the qCDM 
values shown to the LCDM values as Wq changes from —0.5 to —1. Moreover, it is a very 
clean test, because the small-scale S3 is almost independent of all cosmological parameters 
except fi's and Wq. 

To emphasize this point, we show in Fig. |^ S^^s for 3 different qCDM models, with 
different F or erg. They all agree to within a few percent. From eq. (^ and (^, it is not 
hard to see that the normalization of P gets divided out in the combination for S3. However, 
this really refers to the normalization of the nonlinear P. The normalization of the linear P, 
(Tg, should have some effect on S3 through its impact on the shape of the nonlinear power 
spectrum. However, we find that at the small scales which dominate the relevant integrals 
for 5*3, the nonlinear power spectra for most models have rather similar shapes: a slope 
around —1.5 or so. Together with the fact that the nonlinear Q3 is a weak function of scale, 
this explains why the small angular-scale S3 is relatively insensitive to both erg and F. By 
the same reasoning, S3 is quite independent of the spectral tilt as well. Nonlinear clustering 
seems to erase memory of the initial conditions in 5*3, as far as CDM-type power spectra are 
concerned. We therefore have in hand a powerful statistical measure: the small-scale 5*3 is 
almost purely determined by the cosmological energy contents. 



3. Discussion 

How do we understand the cosmological dependence of 6*3? The best way is to go back to 
the definition of the projected mass density k in eq. {^. Observe that w occurs three times 
in the numerator of S3 and four times in the denominator (eq. [Q). This means any overall 
constant multiplying w is going to show up in ^3 (except for Hq/c which is canceled out in 
the combination for k). Among other things, 5*3 scales as l/^m, the reciprocal of the total 



matter content. This means flat matter dominated models generally have lower S3 compared 
to low density models, as is seen in Fig. ^ Moreover, a A or quintessence dominated universe 
has a larger volume out to 2; = 1, compared to an open universe, making w larger and S3 
smaller. Finally, the different fluctuation growth rates in different cosmologies also shift the 
skewness to some extent. We find that the following crude approximation works surprisingly 
well in reproducing our results from integrating eq. @: 

^3- SQs rdxHx)Vrixn9{z)/{l + z)Mx)-'V (10) 









""^ dx[w{xf/r{xfMz)/{l + z)Mxy 



where g{z) is given by eq. (H) or its generalization to include quintessence, Qs is taken to be 
2.7, and n is —1.2 i.e. the power spectrum is assumed to obey a power-law with simply linear 
evolution. The assumptions of a constant Q^ and a power-law P implies a scale-independent 
5*3. The results of applying eq. (p!OD are shown as open circles in Fig. |I|b. 

For more accurate results, we recommend going back to eq. (^, (§) and (|^), which give 
5*3 accurate to within 10% at 6*^ ~ 1' — 5', assuming the sources are at z = 1. We have 
also suggested in §| a useful consistency check: Ak2, a measure that quantifies the degree 
of linearity, should be less than about 30%. For models that are too linear, the hierarchical 
ansatz with a configuration independent Q3 (eq. |]) breaks down. An interesting example is 
provided by the rCDM model simulated by JSW, which is exactly the same as SCDM except 
that r = 0.21 and hence has less power on small scales. As explained before, the highly 
nonlinear 5*3 should be insensitive to F, and our method would predict a rCDM 5*3 very close 
to that of SCDM. JSW found that that the rCDM N-body results are in fact about 30% 
higher than the SCDM results at a few arcminutes. Applying our consistency test shows 
that Ak2 is 2 - 3 times higher for rCDM compared to all other models we have considered. 
Hence, the somewhat large difference between rCDM and SCDM seen in JSW is a reflection 
of their different levels of nonlinearity. Models with as little small-scale power as rCDM are 
probably inconsistent with observations of the Lyman-alpha forest ( [Hui et al. 1997| ; [Croft 
let al. 19981 ). 

We have argued that deep lensing surveys, with a total area of several square degrees 
and background galaxies at 2; ~ 1, should be capable of distinguishing between cosmo log- 
ical models with different energy contents. In particular, contrary to what is indicated by 
perturbation theory, the small-scale skewness can be used to differentiate between curvature 
and cosmological constant dominated models. Moreover, S3 also shows sensitivity to the 
equation of state, Wq, of quintessence models as well, making them distinguishable from A 
models. In practice, however, since we have only one observable in 5*3, one needs to impose 
extra constraints to restrict the range of models when engaging in model testing. For ex- 
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ample, one can assume the class of flat tracker-field models ( [Zlatev et al. 1998 ) where Wq is 



determined by Vtq, and so the only free parameter is Vtq. Another example: making use of 
the fact that for a given fi^ curvature dominated models yield higher S^ than A models, one 
can obtain lower limits on Vim- Additional constraints from other observations such as the 
microwave background and large scale structure are obviously useful. It should also be borne 
in mind that systematic errors, such as those due to the correction of an anisotropic point- 
spread-function, have not been taken into account ( [Kaiser et al. 199"5D . In addition, we have 



assumed that the galaxy redshifts are known, but this is likely achievable by photometric 
techniques. 

A few issues are worth further investigation. A fitting formula for the three-point 
function, which smoothly interpolates between the perturbative and the highly nonlinear 
regimes, could in principle be used to extend our calculation to cover all angular scales. 
At present, no such formula exist for CDM-type spectra (|Scoccimarro fc Frieman 1998| ). 
Moreover, our method can be easily generalized to Sn for arbitrary A^. Such a calculation 
would be useful for the estimation of measurement errors from lensing surveys ( ^coccimarro 



et al. 1999|) . Lastly, our expressions are easily generalizable to measures such as the aperture 



mass ( Schneider et al. 1998| ), which corresponds to using a different smoothing kernel Wq^^ in 
eq. (|^), and its Fourier transform W in the rest of our expressions. Eq. (|1^) should remain 
roughly valid because it is independent of W . 

The author is indebted to Roman Scoccimarro and Matias Zaldarriaga for helpful com- 
ments and for an earlier collaboration which motivated the present work. Special thanks are 
due to Roman Scoccimarro for his very generous help in the course of the investigation. The 
author also thanks Enrique Gaztafiaga for discussions on issues of projection, Martin White 
for discussions on interpretations of N-body simulations, and Zoltan Haiman and Scott Do- 
delson for useful comments. This work was supported by the DOE and the NASA grant 
NAG 5-7092 at Fermilab. 
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Model 


' 'm 


^k 


njnA 


Wq 


F 


0-8 


SCDM 


1 








- 


0.5 


0.6 


OCDM 


0.3 


0.7 





- 


0.21 


0.85 


LCDM 


0.3 





0.7 


-1 


0.21 


0.9 


qCDM 


0.3 





0.7 


-0.5 


0.21 


0.8 



Table 1: A list of models. 
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Fig. 1. — Panel a: A comparison of the skewness obtained from N-body simulations (points 
with error-bars, from 5-10 ray-tracing realizations with no noise added; taken from JSW) 
and from numerical integration of eq. (|^) (lines). The points without error-bars on the far 
left denote, from top to bottom, the perturbation theory predictions of S3 for 0/L/SCDM 
models at Or = 1' (from JSW). Panel b: The lines are the same as before, with the addition of 
predictions for a qCDM model. Points with error-bars denote measurements for 0/L/SCDM 
models from simulated 3° x 3" surveys with 2 x 10^ galaxies per sq. deg. (and random noise 
added; taken from JSW). The open circles on the far left denote the approximate S3 from 
eq. (PUp. All sources/galaxies are assumed to be at z = 1. 
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Fig. 2. — Skewness for three qCDM models. Long-dashed hne - fiducial qCDM as in Fig. ^ 
(Table |ip; dotted line - same qCDM but with F = 0.5; short-dashed line - same qCDM but 
with as = 1.0 



